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ABSTRACT 

Aims. We investigate the origin of the Bry emission of the Herbig Ae star HD 104237 on Astronomical Unit (AU) scales. 



(/} ■ Methods. Using AMBER/VLTI at a spectral resolution R = 1500 we spatially resolve the emission in both the Bry line and the adjacent 
continuum. 

^ . Results. The visibility does not vary between the continuum and the Bry line, even though the line is strongly detected in the spectrum, with a 
peak intensity 35% above the continuum. This demonstrates that the line and continuum emission have similar size scales. We assume that the 
K-band continuum excess originates in a "puffed-up" inner rim of the circumstellar disk, and discuss the likely origin of Bry. 
: Zj Conclusions. We conclude that this emission most likely arises from a compact disk wind, launched from a region 0.2-0.5 AU from the star, 
with a spatial extent similar to that of the near infrared continuum emission region, i.e, very close to the inner rim location. 
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1. Introduction 

The spectra of pre-main sequence stars of all masses show 
prominent strong and broad emission lines, of both hydrogen 
and metals. These lines trace the complex circumstellar envi- 
ronment which characterizes this evolutionary phase, and are 
very likely powered by the associated accretion disks. The 
emission lines are used to infer the physical properties of the 
gas, and to constrain its geometry and dynamics. Their ex- 
act origin, however, is not known. The hydrogen lines, in par- 
ticular, may originate either in the gas which accretes onto 
the s tar from the disk, as in magnetospheric accretion mod- 
els jHartmann et alJfT994l) . or in winds and je ts, driven by the 
interact ion of the accreting disk with a stellar JShu et alTl 994) 
or disk llCasse & Fe rreira 2000) magnetic field. For Herbig Ae 
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stars, it is additionally possible that they form in the inner disk 
dTamhovtseva et al.ll999l). 

For all models, emission in the hydrogen lines is pre- 
dicted to occur over very small spatial scales, a few AUs at 
most. To understand the physical processes that occur at these 
scales, one needs to combine very high spatial resolution with 
enough spectral resolution to resolve the line profile. AMBER, 
the th ree-beam near-IR recombiner of the VLTI ( P etrov etall 
l2003h . simultaneously offers high spatial and high spectral res- 
olution, with the sensitivity required to observe pre-main se- 
quence stars. These capabilities were recently demonstrated by 
Mal bet et alJ J2006I) . who succesfully resolved the luminous 
Herbig Be system MWC297 in both the continuum and the 
Bry line. They showed that the line emission originates from 
an extended wind, while the continuum infrared excess traces 
a dusty accretion disk. 
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Fig. 1. Left: SED predicted by the "puffed-up" rim model used to normalise the K-band continuum visibility. The contributions 
of the A star (long-dashed line), the K3 star (dotted line) and the rim (short dashed line) are sho wn (see the tex t for th e stellar 
parameters). The observations are fr om 2MASS (filled circles), Hipparcos (filled squares), and Malfait et alJ Jl998l) (empty 
circles). An extinction of A v = 0.31 (Ivan den Ancker et alJI 19981) has been used to correct the data for interstellar absorption. 
Center: Comparison of Bry observed with AMBER in the photometric channels (solid line) and ISAAC (dotted line); the dashed 
line shows the ISAAC spectrum smoothed to the spectral resolution of AMBER. Right: Visibility of HD 104237 as a function of 
wavelength. The continuum has been normalised using the star + rim model, as described in the text. The vertical line shows the 
Bry wavelength. 



In this Letter we concentrate on a lower mass, less active 
system, the Herbig Ae system HD 104237. The central emission 
line star, of spectral type between A4V and A8, is surrounded 
by a circumstellar disk, which causes the infrared excess emis- 



sion (Meeus et al 



2001) and drives a jet seen in Ly-a im- 



ages (Gr adv et alJ l2004). The optical spectru m shows a rather 
narro w Ha emission with a P-Cygni profile (Feig elson et alJ 
120031) . The disk is seen almost pole-on (z = 18 Grady 



et al. 200 

-l 



consistent with t he low value o f vsinz (12 km 
s Donati et al. 1 19971) . iDonati et all (1 19971) d etected a stel- 
lar magnetic field of 50 G. Boh m et alJ (2004) have revealed 
the presence of a very close companion of spectral type K3, 
orbiting with a period of ~ 20 days, and which bolometric lu- 
minosity is 10 times fainter than the one of the central star. In 
the near infrared domain, spatially unresolved ISAAC observa- 
tions (Garcia Lopez et al., 2006) show a strong Bry emission 
line, with a peak flux 35% above the continuum. 



2. Observations and data reduction 

AMBER observed HD 104237 on 26 February 2004 on the 
UT2-UT3 baseline of the VLTI, which corresponds to a pro- 
jected length of B = 35m. The instrument was set up to cover 
the [2.121, 2, 21 l]/mi spectral range with a spectral resolution 
of H = 1500, which resolved the profile of the Bry emis- 
sion line at 2.165pm. The data consist of 10 exposures of 500 
frames, with an integration time of 100ms for each frame. The 
integration time is a trade-off between gathering enough flux 
in the photometric channel and preventing excessive contrast 
loss due to fringe motion during the integration. At the time 
of these early observations the star magnitude of K=4.6 was 
close to the sensitivity limit of the medium spectral resolution 
mode, and as a result fringes are only visible in « 25% of the 
frames. In the photometric channels the Bry line is under the 
noise level of the indivual frames, but it is well detected on their 



average and contributes 35% of the total flux (Fig.[0 central 
panel). For comparison we overlay a higher spectral resolution 
spectrum CR - 8900, taken with ISAAC approximately one 
year before the VLTI observations), as well as its smoothing to 
the AMBER spectral resolution. The two spectra are consistent 
with each other, within the combined typical line variability of 
Herbig Ae stars and calibration uncertainties. 

Data reduction followed standard AMBER procedures 

i 1 1 

(Tatulli et al. 2006). To optimize the signal to noise ratio (SNR) 
of the visibility derived from these relatively poor-quality data, 
we pre-selected frames with individual SNR greater than 1.5, 
ensuring that fringes are present in all short exposures that en- 
ter the visibility. 

For these observations the visibility amplitude could not be 
calibrated to an absolute scale using an astronomical calibrator 
(unresolved star), due to non-stationary vibrations in the opti- 
cal train of the UT telescopes. With a 100ms integration time 
these vibrations randomly degrade the fringe contrast of the 
individual frames by a large factor, and the statistics of that at- 
tenuation depends on uncontrolled factors (exact telescope and 
delay line pointing, recent environmental history, etc) that do 
not similarly affect the target star and its calibrator. Fortunately 
the contrast loss from vibrations is achromatic across our small 
relative bandpass, and as a consequence the differential visibil- 
ity, that is, the relative visibility between the spectral channels, 
is unaffected. Our dataset therefore allows to investigate the 
visibility across the Bry emission line compared to that in the 
adjacent continuum. The differential visibilities are accurate to 
approximately 5%. 

The excess near infrared continuum emission in Herbig 

Ae stars most likely originates in the innermost regions of 

their dusty circumstellar disks, at the dust sublimation radius 
i ' — ii 1 i ■ ii 1 



dEisner etalJl2005t llsella etal.j Eo06). To approximately nor- 



malise the visibilities, we thus scaled the observed continuum 
value to the predictions of appropriate theoretical models of 
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the disk inner rim emission dlsella&Nattal l2005h . We com- 
puted the structure o f the "puffed-up" inner rim as described 
in llsella etalJ (E006). We used for the two stars the following 
parameters: T e ff = 8000,^, L* = 3QLq for the primary and 
T e ff = 4700/f, L* = 3Lq for the secondary star, respectively. 
We adopted a distance of D = 1 \5pc, and a disk mass surface 
density of = 2 • 10 3 ■ r 3/2 g/cm 2 (with r in AU). 

Fig.^Oeft) shows that the model using micron size astro- 
nomical silicates produces a very good fit of the SED. The 
star fluxes contribute approximately 30% (respectively 20% 
and 10% for the primary and secondary star, which gives a 
binary flux ratio of ~ 0.5 in the K band) of the total 2.15jum 
flux, and the inner rim appears as a bright ring with radius 
R rim = 0.45AU (3.8mas at D = 115pc). The resulting model 
visibility on the 35m baseline is V — 0.38 (Fig^ right)- Note 
that the orbital period of the spectroscopic binary is 20 days. 
This leads to an average separation of 6s = 0.15AU which 
corresponds to an angular separation of 1.2mas. The central 
system (A star + K3 companion) is therefore completely unre- 
solved on the current baseline, and both stars are located inside 
the dust evaporation radius. 



3. Results and discussions 

Within the fairly small error bars, the visibility does not change 
across the Bry emission line. This result is robust and puts 
strong constraints on the relative spatial extent of the line and 
continuum emission regions, demonstrating that they have very 
similar apparent sizes. We use this constraint to probe the pro- 
cesses responsible for the the Bry emission in this star, and con- 
sider in turn the three main mechanisms usually invoked to in- 
terpret the hydrogen line emission in pre-main-sequence stars. 
We translate each mechanism to a simple geometrical models 
of specific spatial extension, with the line strength fixed at the 
observational value, and evaluate the resulting visibility across 
the line 1 . 

Magneto spheric accretion: in such a model matter in- 
falls on the star along magnetic field lines, and the base of 
this i nfalling flow is (appro ximately) inside the corotation ra- 
dius jMuzerolle et alJl2004 . F or HP 104237 we fin d R corol = 
0.07AU, using vsin; = 12km/slDonati et alJ J 19971) and an in- 
clination of i — 18 JCradv et all2 004). Note that this value is, 
at our spatial resolution level, similar to the separation of the 



1 Note that this analysis is valid as long as the continuum emis- 
sion is resolved. In our case the continuum is calibrated by a model 
and not by an unresolved star, and therefore there might be a chance 
that the visibility of the continuum is close to 1. However, this pe- 
culiar scenario appears very unlikey since the puffed-up rim model 
has been shown to be well representative of the very close environ- 
ment of Herbig Ae stars over a large luminosity range iMonnier et alJ 
120051) . Furthermore, the hypothesis of a resolved continuum emission 
for HP 10 4237 is stro ngly supported by its measured accreti on rate of 
10- 8 M o iGradv et alJl2004 . Indeed. iMuzerolle et alJ J200I showed 
that for weak accretion rates (M < lCT 7 Mo/yr), the gaz accreting onto 
the star is optically thin. It therefore does not shield the dust from the 
stellar radiation and consequently, the star must exhibit an inner re- 
gion free of dust, large enough to be resolved by the interferometer at 
the 35m baseline. 
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Fig. 2. Comparison between the observed visibilities (empty 
square with error bars) and the predictions (solid curves) of the 
simple geometrical models for the Bry emission (sketched in 
the same panels). The observed visibilities are scaled to match 
the continuum value predicted by the "puffed-up" inner rim 
model (Isella & Natta, 2005) as described in Sec. 2. The contin- 
uum emission arises both from the stellar photosphere («20%) 
and from the dusty disk inner rim, located at the dust evapora- 
tion distance R rml = 0.45AU and which appears as the bright 
gray scale ring. The Bry emission regions are shown as grid 
surfaces. The three panels illustrate the different models dis- 
cussed in Sec. 3: the upper panel represents the magneto spheric 
accretion model in which the Bry emission originates very 
close to the star, inside the corotational radius R corot =0.07AU; 
in the middle panel the Bry emission originates between R corot 
and the rim radius R rm =0.45AU, representing the gas within 
the disk model; the bottom panel shows the outflowing wind 
model, in which the emission is confined close to the inner rim, 
between ~ 0.2AU and ~ 0.5AU. 



binary. Adopting R slar and R coro , as the inner and outer limits 
of the magnetospheric accreting region, Fig. |2 (upper panel) 
demonstrates that Bry emission is then confined much closer 
to the star than the dusty rim. The predicted visibility therefore 
increases significantly in the line, contrary to the observational 
result. Explaining the observed visibility with magnetospheric 
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accretion requires the corotation radius to approach the inner 
rim, which would need an unrealistically lower stellar rota- 
tional velocity (v < 2km/s). The magn etic field of HD 104237 
is weak, B — 5QG ( Do nati et al Jl 19971) . and the corresponding 
magnetospheric truncation radius (R mag = 0.01 8 AU, Shu et al. 

is smaller than the corotation radius. Replacing R CO rot 
by Rmag as the outer radius of the magnetospheric accreting re- 
gion would thus only reinforce our conclusion. Magnetospheric 
accretion therefore cannot be responsible for most of the Bry 
emission. 

Gas in the disk: sinc e HP 104237 accretes matter (a few 
10" 8 M Q /yr, Grady et al. l2004 Garcia Lopez et al. 200§, an 
optically thin gaseous disk (Muzerolle et al. 2004) must extend 
inward from the dust evaporation rim to either the magneto- 
spheric truncation radius or the corotation radius; this disk re- 
gion is ionised and may, in principle, emit the observed Bry 
line. To model this scenario, we assume a constant Bry sur- 
face brightness between R coror (or R mag ) and R r t m . The resulting 
visibilities (Fig.fJl middle panel) are also inconsistent with the 
observations, though not quite as severely. 

Outflowing wind: the remaining possibility is that Bry is 
emitted in a wind or jet. The bulk of the emission is then con- 
fined to the regions of highest gas density, i.e., to the base of 
the wind/jet. Since the wind/jet must be seen almost pole-on, 
we assume in our toy model that emission is confined to a 
ring of width AR, with AR/R/ ~ 0.5, where R, is the ring in- 
ner radius, that is the wind launc hing point. This as sumption 
is guided by the wind models of Na tta et alJ Q988I) . and by 
the more r ecent simulations of the Pa/3 emission in jets of T 
Tauri stars jThiebaut et alJ2003h . For pole-on outflows most of 
the intensity originates in a ring with these approximate prop- 
erties. Fixing AR/Ri = 0.5 and adjusting /?,-, this toy model 
reproduces the interferometric data for /?, between ~ 0.25 AU 
and 0.35 AU (Fig. |3 bottom panel). Allowing different rela- 
tive widths (AR from 10% to 100% of R { ), we find that the Bry 
emission has to originate between ~ 0.2 and 0.5 AU. Our in- 
terferometric measurements are therefore consistent with the 
Bry emission coming from the base of a wind, originating in a 
disk region close to the location of the dusty rim. The presence 
of a wind in HD 104237 is well establishe d on (much) l arger 
scales, by the Lya bipolar microjet (Gradvetal. 2004) and 
by the P-Cygni profile of the Ha line (Feigels on et alJ l2003). 
Given our oversimplified geometrical models, it would be in- 
appropriate to speculate too much on the detailed nature of the 
wind. We note however that the launching region of an X-wind, 
driven by the stellar magnetic field (Shu et al. 1994), is close 
to the corotation radius and too small to be consistent with the 
present data. Making the X-wind acceptable needs relaxing the 
assumption that most of the line emission originates near the 
launching point, and instead having the brightest Bry region a 
factor 5-8 further away from the star. The disk-wind scenario 
by contrast has its launching point a few tenths of AU from the 
star and needs no tweaking. This preference for a disk-wind is 
consistent with the size of jet launching reg ions inferred from 
the rotation of TTauri jets JCoffev et al.l2004l) . 

It is important to recall here that our analysis assumes that 
the continuum originates in a dusty "puffed-up" ring. That 
model has been shown to match Herbig Ae and T Tauri stars 



over a large luminosity range JMonnier et alJEo05i) . but it ob- 
viously needs to be verified for the specific case of HD 104237, 
by obtaining calibrated interferometric observations over a few 
baselines. 

4. Conclusions 

We have presented interferometric observations of the Herbig 
Ae star HD 104237, obtained with the AMBER/VLTI instru- 
ment with H = 1500 high spectral resolution. The visibility 
is identical in the Bry line and in the continuum, even though 
the line represents 35% of the total 2.165yum flux. This imme- 
diately implies that the line and continuum emission regions 
have the same apparent size. 

Scaling the continuum visibility with a "puffed-up" inner 
rim model, and using simple toy models to describe the Bry 
emission, we have shown that the line emission is unlikely to 
originate in either magnetospheric accreting columns of gas or 
in the gaseous disk. It is much more likely to come from a com- 
pact outflowing disk wind launched in the vicinity of the rim, 
about 0.5 AU from the star. This does not preclude accretion 
from nonethel ess occuring al ong the stellar magnetic field de- 
tected bv lDonati et alJ Jl997l) . and accreting matter might even 
dominate the optical hydrogen line emission, but our observa- 
tions show that the bulk of the Bry emission in HD 104237 is 
unlikely to originate in magnetospheric accreting matter. 

Our results show that AMBER/VLTI is a powerful diagnos- 
tic of the origin of the line emission in young stellar objects. 
Observations of a consistent sample of objects will strongly 
constrain the wind launching mechanism. 
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